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Special solutions of the LTB family representing collapsing over-dense regions extending asymp- 
totically to an expanding closed, open, or flat FRW model are found. These solutions may be con- 
sidered as representing dynamical mass condensations leading to black holes immersed in a FRW 
universe. We study the dynamics of the collapsing region, and its density profile. The question of 
the strength of the central singularity and its nakedness, as well as the existence of an apparent 
horizon and an event horizon is dealt with in detail, shedding light to the notion of cosmological 
black holes. Differences to the Schwarzschild black hole are addressed. 

PACS numbers: 95.30.Sf,98.80.-k, 98.62. Js, 98.65.-r 



I. INTRODUCTION 

Let us use the term cosmological black hole for any 
solution of Einstein equations representing a collapsing 
overdensity region in a cosmological background leading 
to an infinite density at its center There have been 
different attempts to construct solutions of Einstein 
equations representing such a collapsing central mass. 
Gluing of a Schwarzschild manifold to an expanding 
FRW manifold is one of the first attempts to construct 
such a cosmological black hole, as done first by Einstein 
and Straus 0. Depending of the way the model is 
constructed, jane is led to un-physical behavior of the 
trajectories [3j]. 

Models not based on a cut-and-paste technology 
is much more interesting giving more information on 
the behavior of the mass condensation within a FRW 
universe model. The first attempt is due to McVittie 
introducing a spacetime metric that represents a 
point mass embedded in a Friedmann-Robertson- Walker 
(FRW) universe. There have been many other attempts 
to construct cosmological black holes such as Nolan 
interior solution [1], and Sultana-Dyer solution[l|, each 
of them contrasting some of the features one expect 
from theory or observation. 

The interest for cosmological black holes in the past 
has been mainly from the theoretical side to understand 
concepts like black hole, siiigularity, horizon, and ther- 
modynamics of black holes [q. Indeed, the conventional 
definition of black holes implies an asymptotically flat 
space-time and a global definition of the event horizon. 
In practice, however, the universe is not asymptotically 
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flat. The need for local definition of black holes and 
their horizons has led to concepts such as Hayward's 
trapping horizon Q, Ashtekar's isolated horizon 
Ashtekar and Krishnan's dynamical horizon llOll . and 
Booth and Fairhurst's slowly evolving horizon 
There are cases where both apparent and event horizon 
maybe possibly defined. For example, for dynamical 
black holes one may define the event horizon as the very 
last ray to reach future null infinity or the light ray that 
divides those observer who cannot escape the future 
singularity from those that can [2^. Eardley proposed 
the conjecture that in such cases trapped surfaces can 
be deformed to get arbitrarily close to the event horizon 
p^ . Numerical evidence was provided in [ip] and later 
proved analytically for the Vaidya metric |14| . 

The precision cosmology has opened a new arena 
for questions like cosmological black holes and their 
behavior. New observation of our galactic center allow 
to resolve phenomena near the Schwarzschild horizon 
of the central black hole 0] • It is therefore desirable to 
have black hole models embedded in cosmological envi- 
ronment to see if there may be considerable differences 
to the familiar Schwarzschild black hole. There have 
been also increasing interest in the gravitational lensing 
by a cosmological mass condensation such as a cluster 
of galaxies in a cosmological background. The simplest 
cases are the Kottler and the Einstein-Straus model 
[2^ . The more complex situation is lensing by a mass 
condensation within a dynamical background. 

Now, a widely used metric to describe the gravita- 
tional collapse of a spherically symmetric dust cloud 
is the so-called Tolman-Bondi-Lemaitre (LTB) metric 
[13] ■ These models have been extensively studied for 
the validity of the cosmic censorship conjecture [l8|,[l9j 
and [lO]. In particular, we know already [21[ that, 
depending upon the initial conditions defined in terms 
of the initial density and velocity profiles from which the 
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collapse develops, the central shell-focusing singularity 
at r = can be either a black hole or a locally or 
globally naked singularity. We may note however, that 
in all these papers a compact LTB region is glued to 
the Schwarzschild metric or the FRW outer universe 
[2^ . Therefore, the results have to be taken cautiously: 
any principally existent event horizon is cut off by the 
outer static or homogeneous space-time. The statement 
may still be correct that in a dynamic spacetime the 
cosmic censorship hypotheses is valid, as discussed in 
[2^ . It is also possible to glue two different LTB metrics 
to study the structure formation out of an initial mass 
condensation or the formation of a galaxy with a central 
black hole and [11]. Here again the structure of 
the metric outside the mass condensation is defined by 
hand to match with a specific galaxy or cluster feature. 
Faraoni et al have tried to change McVitte metric so 
that it resemble a collapsing mass condensation. Their 
solution, however, represents a singularity within a 
horizon embedded in a universe filled with a non-perfect 
fluid where the change of the mass is not because of 
the in- falling matter but the heat flow [23|. This metric 
gives us no clue whatsoever about the dynamics of a 
possible collapsed mass condensation. Harada et al, 
being interested in the behavior of primordial black holes 
within cosmological models with a varying gravitational 
constant, use a LTB solution to study the evolution of 
a background scalar field when a black hole forms from 
the collapse of dust in a flat Friedmann universe probing 
the gravitational memory [Tsj . 

Our goal is to look for a model of a cosmological black 
hole, i.e. a mass condensation leading to a singularity 
within a FRW universe universe. In this paper we 
propose the models for closed-, open-, and flat FRW 
universe studying their density profiles, singularities, 
and horizon behaviors. There are many nontrivial ques- 
tions to be answered before understanding in detail the 
differences of these cosmological black holes to the famil- 
iar Schwarzschild ones, which are beyond the scope of 
this paper and are to be dealt with in future publications. 

The question of singularities and the definition of a 
black hole in such a dynamical environment has been sub- 
ject of different studies in the last 15 years. We review 
very shortly different definitions of horizons in section II 
as a reference to the properties of model solutions we pro- 
pose. Some initial attempts to model black holes within 
a FRW universe is introduced in section III. Section IV is 
devoted to the LTB metric as the generic solution repre- 
senting a spherically symmetric ideal fluid. Section V is 
devoted to different models of cosmological black holes, 
their dynamics, density profile, apparent and event hori- 
zons, and singularities. The question of strength and the 
nakedness of singularities are dealt with in section VI. 
We then conclude in section VII. Throughout the paper 
we assume SnG = c = 1. 



II. LOCAL DEFINITIONS OF BLACK HOLES 

Standard definition of black holes Q needs some 
global assumptions such as regular predictability and 
asymptotic flatness. In the cosmological context con- 
cepts of asymptotic flatness and regular predictability 
have no application. This has already been noticed by 
Demianski and Lasota [s^ stressing the fact that in 
the cosmological context the standard global definition 
of black holes using event horizons may not be used 
any more. Tipler [3J] also present a definition of black 
hole in non asymptotically flat space time, but these 
definition did not have comprehensive property of black 
hole such as thermodynamic laws. In the last decade 
the interest in a local definition of black holes has led to 
four different concepts based primarily on the concept 
of the apparent horizon. 

Let us start by assuming a spacelike two surface 
S with two normal null vectors and n° on it. The 
corresponding expansions are then defined as ()[g),0[n)- 

Definition 18]. A trapping horizon H is a hy- 
persurface in a 4-dimensional spacetime that is foliated 
by 2-surfaces such that 6'(£) \h= 0, 6'(„) |//^ 0, and 
£nQ{i) 0. A trapping horizon is called outer if 

£nd{i) \h< 0, inner if £„6'(£) \h> 0, future if 6'(„) \h< 0, 
and past if 6'(„) \h> 0. The most relevant case in the 
context of black holes is the future outer trapping hori- 
2on(FOTH). 

Definition 2 Q. A weakly isolated horizon is a three- 
surface H such that : 

1. H is null; 

2. The expansion 6'(£) where being null and 
normal to the foliations S oi H; 

3. — Tq^" is future directed and causal; 

4. £iUJa ~ 0, where usa — —ni,Wai^, and the arrow 
indicates a pull-back to H. 

Weakly isolated horizon is a useful term to be used for 
characterization of black holes not interacting with their 
surroundings, and corresponds to isolated equilibrium 
states in thermodynamics. These definition do not ap- 
ply to cosmological mass condensations because of their 
dynamical behavior. 

Definitions lOj. A marginally trapped tube T (MTT) 
is a hypersurface in a 4-dimensional spacetime that is 
foliated by two-surfaces S, called marginally trapped 
surfaces, such that 0(„)|t < and ^(£)|t = 0. MTTs 
have no restriction on their signature, which is allowed 
to vary over the hypersurface. This is a generalization 
of the familiar concept of the apparent horizon p^ . 
If a MTT is everywhere spacelike it is referred to as 
a dynamical horizon. If it is everywhere timelike it 
is called a timelike membrane (TLM). In case it is 
everywhere null and non-expanding then we have an 
isolated horizon. The apparent horizons evolving in the 
our proposed models will not be everywhere spacelike 
and will have a complex behavior. 
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Irrespective of different concepts related to the ap- 
parent horizon we may still compromise on a defini- 
tion of event horizon differing principally from the ap- 
parent horizon. We follow the definition of [1^ as the 
very last ray to reach future null infinity or the light 
ray that divides those observer who cannot escape the 
future singularity from those that can. We will see in 
the next sections that cosmological black holes may have 
distinct apparent and even horizons, in contrast to the 
Schwarzschild black hole. 



III. EXISTING METRICS REPRESENTING 
OVER-DENSITIES WITHIN A COSMOLOGICAL 
BACKGROUND AND THEIR DEFICIENCIES 

A. McVitties solutions 



black hole is primordial in the sense that it forms ab ini- 
tio with the big bang singularity and therefore does not 
represent the gravitational collapse of a matter distribu- 
tion. 

This metric is given by 

ds^ = t%\ )dt^ dtdr-{\ )dr^ -r^d?L\ 

r r r 

Though the metric has the same causal characteristics 
as the Schwarzschild spacetime, there are significant dif- 
ferences for timelike geodesies. In particular an increase 
in the perihelion precession and the non-existence of cir- 
cular timelike orbits should be mentioned. The matter 
content is described by a non-comoving two-fluid source, 
one of which is a dust and the other is a null fluid. At 
late times the dust becomes superluminal near horizon 
violating the energy condition. 



In 1933, McVittie found an exact solution of Ein- 
steins equations for a perfect fluid mimicking a black hole 
embedded in a cosmological background. McVitties so- 
lutions can be written in the form 



ds' 



1 _ M. 

■'^ fde+eP^'\l+—Y{dr^+h^d^f ), (1) 



M 



2N 



2N' 



where M = me'^'*^/^ and m is a constant. Functions 
h(r) and N(r) depend on a constant fc, and are given, 
respectively, by 

sinh(r) k = 1 
h{r) = I r fc = 

sin(r) fc = — 1 

2sinh(r/2) fc 1 
N{r) ^ { r k = (2) 

2sin(r/2) fc = -1 

This metric represents a point mass embedded into an 
isotropic universe. It possesses a curvature singularity at 
proper radius R — 2m, in contrast to the Schwarzschild 
metric where there is a coordinate singularity. It has been 
shown that this singularity is space- like and weak[35j. 
The interpretation of the metric in the region R < 2m 
is also not clear |35i]. Therefore, the McVittie's metric is 
not a suitable solution of Einstein equations to represent 
the collapse of a spherical mass distribution with over- 
density within a cosmological setting. 



B. Sultana-Dyer solution 

Recently Sultana and Dyer [l| found an exact solu- 
tion representing a primordial cosmological black hole. 
It describes an expanding event horizon in the asymp- 
totic background of the Einstein-de Sitter universe. The 



IV. INTRODUCING REALISTIC MODELS OF 
COSMOLOGICAL MASS CONDENSATION 

There maybe different ways of constructing solutions 
of Einstein equations representing a collapsing mass con- 
centration in a FRW background, as the preceding sec- 
tions show. We choose the direct way of a cosmological 
spherical symmetric isotropic solution, and look for an 
overdensity mass distribution within the model universe 
undergoing a collapse to see if and how a singularity rep- 
resenting a black hole emerges. To begin with, we choose 
a so-called flat LTB metric. This is the simplest spheri- 
cally symmetric solution of Einstein equations represent- 
ing an inhomogeneous dust distribution [T7| . 



A. LTB metric 

The LTB metric may be written in synchronous coor- 
dinates as 



ds^ = dt^ 



R'- 



-dr^ - R{t,r)^dn^ 



1 + fir) 

It represents a pressure-less perfect fluid satisfying 



2M'{r) 
R^R' 



R^ = .f + 



2M 
~R' 



(4) 



(5) 



Here dot and prime denote partial derivatives with re- 
spect to the parameters t and r respectively. The angular 
distance i?, depending on the value of /, is given by 

M 

R = — ^(1 - cos77(r, t)), 
77-sinr;= ^^j^(<-t„(r)), (6) 



R = (-/) 



1/2 sin{ri) 
1 — cosrj ' 



(7) 
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for / < 0, and 

i?=(^Af)^(i-i„)§, (8) 

for / = 0, and 

M 

R = — (cosh7y(r, t) - 1), 

J3/2 

sinh77-r;== -^(t-t„(r)), (9) 

for / > 0. 

The metric is covariant under the rescahng r r{r). 
Therefore, one can fix one of the three free parameter 
of the metric, i.e. tn{r), f{r), and M{r). The function 
M(r) corresponds to the Misner-Sharp mass in general 
relativity, as shown in the general case of spherically sym- 
metric solutions of Einstein equations [l6j. 
There are two generic singularities of this metric: the 
shell focusing singularity at R{t, r) = 0, and the shell 
crossing one at R'{t,r) = 0. To get rid of the complex- 
ity of the shell focusing singularity, corresponding to a 
non-simultaneous big bang singularity, we will assume 
tn{r) — 0. This will enable us to concentrate on the be- 
havior of the collapse of an overdensity region in an ex- 
panding universe without interfering with the complexity 
of the inherent bang singularity of the metric. 
Now, an expanding universe means generally i? > 0. 
However, in a region around the center it may hap- 
pen that i? < 0, corresponding to the collapsing region. 
It is then easy to show that in this collapsing region 

/ 2M I f' / 2M I f 

« (1 - ^vfrr ^(") ^ " ^7^) < ^- 

fore, R = 2M, is obviously a marginally trapped tube, as 
defined in section 2, representing an apparent horizon ac- 
cording to the familiar definitions [gI. [io| . It will turn out 
that this apparent horizon is not always spacelike and 
can have a complicated behavior for different r, as was 
first seen in [36| . 

B. Behavior of the curvature function /(r) 

Now, we are interested in an expanding universe, 
meaning generally R > 0. However, in a region around 
the center we expect to have a late time behavior R < 
corresponding to the collapse phase of the overdensity 
region. From equations ([6]), ([8]), and ([9]) we infer that 
to have a collapsing region one has to ask for f{r) < 
in that region. In contrast, the universe outside the 
collapsing region being expanding leave us to choose 
/(r) > 0, /(r) = 0, or /(r) < depending on the model. 
We may have an asymptotically flat FRW universe, 
however, with /(r) > or /(r) < tending to zero for 
large r. 

Now, we have still to make a choice for /(r) at the 
center r = 0. Expecting the mass M to be zero at ?■ = 
to avoid a central singularity, we see from ^ and ^ 



that ^jj-\r=o — const, or /(r = 0) = [s^]. This give us 
different possibilities of the function /(r) to behave as 
shown in Fig[TJ 



f(r) 




FIG. 1: Different behaviors of the curvature function /(r). 



V. CONSTRUCTION OF MODELS 

We have now the necessary prerequisites to construct 
our models of mass condensation immersed in FRW mod- 
els leading to singularities and representing cosmologi- 
cal black holes. Our cosmological black hole solutions 
evolve from mass condensations within closed, open, or 
flat FRW universes, leading to singularities having differ- 
ent horizons, and providing us examples of collapsed re- 
gions behaving differently to known Schwarzschild ones. 

A. Example I: / < 0: asymptotically closed LTB 
metric 

As mentioned before, we are free to choose one of the 
three parameters of the LTB metric. Assuming a nega- 
tive f{r), we may choose r such that /(r) — —M{r)/r 
(Tst . Now, let us choose the mass function AI such that 

^ ' l + r^ 

where a and a are constants to be defined properly. We 
then obtain from ^ 

R = r(l — cosrj{r, t)) 

77 - sin[r^) = ^m±^2^/^at. (10) 

We are free to fix a and a such that for the present 
time, to I the region around the center of the overdensity, 
r = 0, is collapsing while far from the center the universe 
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expands. Note that in contrast to the famiUar FRW uni- 
verse, where the scale factor as a function of time, t, is 
an exphcit function having a straightforward behavior. 
In the LTB case, R{t) playing the role of the scale factor 
is an implicit function of time and comoving coordinate 
r given by ([6]). We now fix a and a such that r = 
corresponds to 77 = ^ , and r ^ 1 corresponds to = 
(Fig. [2]) . We then find a ~ 2^15 with to being the present 
time, and a ~ 7. 

Now, the expansion phase of the model is given hy R 
We then see from ([7]) that the region around r = 0, cor- 
responding to 77 ~ is always collapsing for any time 
t, while the regions far from the center, r >> 0, at the 
present time, corresponding to ~ are expanding. 
Note that this bound LTB model, similar to the closed 
FRW one, has a maximum comoving radius correspond- 
ing to /(r) = -1. 




Big Crunch singularity 



FIG. 2; Evolution of the Cauchy surfaces. 

The density evolution and the causal structure of the 
model is shown in Fig. [3l We see clearly how the central 
overdensity region collapses to a singularity at ?■ = 0, 
while the universe is expanding. Note also how the slope 
of outgoing null geodesies tend to infinity in the vicinity 
of the singularity, i.e. R' +00 at R = 0. 
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FIG. 3: The case of the asymptotically closed universe: in 
the central region the density increases with time indefinitely 
while far from the center the density is decreasing with time. 
The apparent horizon and the trapped region is shown in the 
lower diagram. 



where a is a constant having the dimension [a] = 
We fix a by ato = 37r/2. Similar to our previous model 
I, this value of a corresponds to the collapsing mass 
condensation around r = starting in the expanding 
phase of the bound LTB model. 



B. Example II: / < 0, lim,.^oo /(r) -» 0; 
asymptotically flat LTB metric 1 

Our favorite choice is a solution representing a collaps- 
ing overdensity region at the center and a flat FRW far 
from the overdensity region. Of course the overdensity 
region may take part in the expansion of the universe 
at early times but gradually reversing the expansion and 
start collapsing. To achieve this, we require /(r) < 
and /(r) — > when r ^ 00. This choice give us trivially 
Af (0) = 0. 

Let us now make the ansatz j(r) = — re^'' leading to 

M(r) = irV2(l + r3/2), 
a 



Equation ([6]), ([8]) then leads to 

R^ ^ ^ \ ' {l-cos7j{r,t)), 
ae ' 

We have plotted the density evolution and casual struc- 
ture of this model in Figd) 

As a result of R' +00 near the singularity, the slope 
of the outgoing null geodesies becomes infinite at the cen- 
tral singularity. Again we see clearly how the collapse of 
the central region and the evolution of the apparent hori- 
zon separates the overdense region from the expanding 
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'Singularity 
Apparent horizon 



Trapped region 



Event horizon 



Untrapped 
region 



FIG. 4: The density profile for the cosmic blaclc hole within a 
closed but asymptotically flat universe. The causal structure 
is shown below. Note the behavior of the event horizon for 
arbitrary large but finite t. 



universe. 

The negativity of the curvature function f(r) means 
that, although the universe is asymptotically flat, wait- 
ing enough, every slice r = constant will collapse to the 
central region. We may , however, define an event hori- 
zon according to the definition of section 2 for any large 
but finite time, as shown in Figd) 



C. Example III: / > 0, /(r) when r oo; 
asymptotically flat LTB metric 2 

What would happen if we choose the curvature func- 
tion /(?■) such that it tends to zero for large r while it 
is positive? We still have a model which tends to a flat 
FRW at large distances from the center, but having a 
density less than the critical one. 

Let us make the ansatz /(r) — —r{e~^— ^J, with n = 2 
and c — 20000, leading to 



M{r) = -rV2(i + r3/2)^ 



where a is a constant having the dimension [a] — [L]^"^ . 
We fix a by requiring aio = 37r/2. Equation ©, ^ then 
leads to 



r''+20000^ 



77 — smrj 



and for / > region. 



r-^+20000 



1.5 



-at. 



V^(l + r3/2) 



(l + r3/2) 



(cosh77(r, t) — 1), 



(12) 



-1-20000 



rj — sinhrj = 



r^+20000 



-e-n 



-r^l.5 



(l-^r3/2) 



-at. 



(13) 



The solution is continuous at r = 1, as can be checked 
by evaluating R, R' , R' , R, and R' a.t r = l(see the 
appendix) . 

We have plotted the density evolution and casual struc- 
ture of this model in Fig[5] 

The term -prf^ is responsible for f{r) being positive 
and tending to zero for large r given n > 2 and c >> 1. 
Let us check if this may cause shell crossing in the region 
where f'{r) < while / > 0. Using ^ we obtain 



R' _M' 

'r ~ 



f 3 



(14) 



where 



< 



/ sinhrj {sinhrj - 

"4 — (coshrj-iy 



< 1. The condition for 
Min < For 



no shell crossing singularity is then ^ :i04-i 

04 ~ 1, corresponding to77>> 1 ort>> 1 the inequal- 
ity breaks down leading to a shell crossing singularity. 
The shell crossing, however, can be shifted to arbitrary 
large t b y ch oosing /' << 1 corresponding to n >> 1 and 
c >> 1 [33] ■ Therefore, for the model we are proposing 
the shell crossing will happen out of the range of appli- 
cability of it. 

As a result of R' +00 near the singularity, the slope 
of the outgoing null geodesies become infinite at the cen- 
tral singularity. Again we see clearly how the collapse 
of the central region and the evolution of the apparent 
horizon separates the overdense region from the expand- 
ing universe. There is an event horizon defined by the 
very last ray to reach future null infinity and separates 
those observer who can not scape the future singularity 
from those that can. A fixed r — ro value, being the non- 
trivial root of /(r) = 0, divides the absolute collapsing 
region from the absolute expanding region. We may be 
living in a region inside the event horizon but outside the 
apparent one without noticing it soon! 

This solution represents a collapsing mass within an 
asymptotically flat FRW universe. The collapsed region 
is dynamical in the sense that its mass is not constant. 
In fact the rate of change of the Misner-Sharp energy is 

given by ^^^\R=const = ^^^^^\B^=const > because 
> 0, R'dr + Rdt = 0, i?' > 0, and i? < for 
collapsing region, so ^\R=const > 0. Therefore, it is 
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Untrapped 
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Singularity 
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FIG. 5: Evolution of the cosmic black hole within an open but 
asymptotically flat universe is similar to the closed case. The 
causal structure, however, is significantly different, as seen 
from the lower diagram. Result of the numerical calculation 
of the locations of the event horizon, apparent horizon and 
the singularity is also shown. 



D. Example IV: / > 0: asymptotically open FRW 
metric 

Now we look for a solution which goes to an open FRW 
metric at distances far from the center. At the same time 
one should take care of the conditions M{0) = and 

M{0) 



7^ oo. Let us choose 
fir) = - 



and 



M(r) = ir3/2(l + r3/2), 



where a is a constant, which may be fixed by assuming 
r = at the present time to corresponding to 77 = 
This leads to ato = 37r/2 + 1 . We then obtain from ^ 

V^(l + r3/2) 



R = 



a{l - r) 



-(1 - C0S77(r, t)), 



> , (l-r)3/2 
77 — sin[r\) = — TTH^at, 



for r < 1, and 



R 



^3/2 



a{r - 1) 



(l + r3/2) 



{coshr]{r, t) — 1) 



(15) 



(r - 1)3/2 
sinhrj — 77 = — TTTTzat, 



1 



,3/2) 



(16) 



for r > 1. The solution is again continuous at r = 1, 
as can be checked by evaluating R, R' , R' , R, and R' at 
r — l(see the appendix). 

The resulting density profile and the causal structure 
is plotted in FiglH Obviously a singularity at the ori- 
gin forms gradually while the universe is expanding. The 
causal structure is also similar to the open but asymp- 
totically flat case. 

This solution represents a collapsing mass within 
an open FRW universe. The collapsed region is again 
dynamical in the sense that its mass is not constant, and 
the rate of change of the Misner-Sharp energy is given 
by the same amount as the previous model. Therefore, 
concepts of isolated horizon and slowly evolving horizon 
do not apply to this case. 



VI. CHARACTERISTICS OF SINGULARITIES 
OF PROPOSED MODELS 



clear that concepts such as isolated horizon and slowly 
evolving horizon do not apply to this case. 



We have avoided in the models proposed the shell 
crossing singularities except example III with a late time 
shell crossing singularity. 

The shell focusing singularities, however, are un- 
avoidable and in fact it is what we are looking for to 
study characteristics of cosmological black holes. An 
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Singularity 

► Apparent horizon 



Trapped region 
Event liorizon 




FIG. 6: The case of asymptotically open model: the density 
profile is similar to the previous open but asymptotically flat 
case, except for the less mass concentrated in the central re- 
gion. The causal structure is also similar. The locations of 
the event horizon, apparent horizon and the singularity are 
also calculated numerically. The separation between the sin- 
gularity and the apparent horizon is not clear here due to the 
scale chosen. 



important asp ect of such a singularity is its gravitational 
strength [31[ . which is an important differentiating 
feature of black hole. 



A. Strength of the shell focusing singularities 

Heuristically, a singularity is termed gravitationally 
strong, or simply strong, if it destroys by crushing or 
stretching any object which falls into it. The prototype 
of such a singularity is the Schwarzschild one: a radially 
infalling object is infinitely stretched in the radial direc- 
tion and crushed in the tangential directions, with the 
net result of crushing to zero volume. Otherwise a sin- 
gularity is termed weak where no object falling into it is 
destroyed. To check the strength of singularities of our 



31| 



models we use the criteria defined by Clarke 
Let k'^ be the tangent vector to the ingoing null geodesic, 
and A the corresponding affine parameter being zero at 
the center. R^^i, being the Ricci tensor, the singularity is 
said to be strong if 



(17) 



For a 
kt^k^R 



general 



LTB 

M' 



metric one obtains easily 
For the three interesting cases 



of cosmological black holes in flat and open LTB models 
we have done the calculation along the lines of the 
using appropriate coordinates near singularity. For our 
cases (lllj|13IIT5)) we obtain after some calculation = 
for r — > 0. Therefore, shell focusing singularities occur- 
ring in the center of the models we are proposing are 
week. This is in contrast to the Schwarzschild singular- 
ity which is a strong one. We leave it to future studies if 
this weakness is generic of any cosmological black holes. 



B. Nakedness of singularities 

We know already from Oppenheimer-Snyder collapse 
of a homogeneous dust distribution how the shells be- 
come singular at the same time, and thus none of them 
crosses. In the case of spherically symmetric inhomoge- 
neous matter configurations, however, the proper time 
of collapse depends on the comoving radius r. Thus the 
piling up of neighboring matter shells at finite proper ra- 
dius can occur, thereby producing two-dimensional caus- 
tics where the energy density and some curvature compo- 
nents diverge. These singularities can be locally naked, 
but they are gravitationally weak [3, HI], i.e. curva- 
ture invariants and tidal forces remain finite. It has also 
been shown that analytic continuations of the metric, in 
a distributional sense, can always be found in the neigh- 
borhood of the singularity . 

Models proposed in this paper are, however, free from 
shell crossing singularities. The shell crossing singularity 
of example III at late times does not influence the follow- 
ing argumentation. Conditions for the absence of shell 
crossing singularities have been studied in detail in [STf . 
In our case these conditions are equivalent to M'(r) > 
and R' > 0, which are satisfied by the models discussed 
above. We may then conclude that 
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dt I 
dt I 



= 1- 



2M' 



< 1. 



(18) 



Therefore, the condition for the apparent horizon R = 

— I d 

2M to be spacehke is, i.e. -1 < ^f,^ < 1, leads to 

the condition R' — M' > 0, which is not everywhere sat- 
isfied in our modeL As a result we notice that apparent 
horizons of the models proposed here are not spacelike 
everywhere. Such a behavior has already been discussed 
in[36|. 

The case of shell focusing singularities is, however, a dif- 
ferent one. Irrespective of the behavior of the apparent 
and event horizons, it is then a relevant question if the 
shell focusing singularity could be a naked one. We no- 
tice that the slope of the outgoing null geodesies at the 
singularity are greater than the slope of the singularity 
itself. Therefore, the singularity is spacelike and no time- 
like or null geodesic can come out of the singularity. We 
then conclude that the singularities we are facing can not 
be naked. 



VII. DISCUSSION AND CONCLUSIONS 

Unlike models discussed so far in the literature, we 
have constructed models of mass condensation within the 
FRW universe leading to cosmological black holes with- 
out having the usual pathologies we know from other 
models: the cosmic fluid is dust and ideal producing a sin- 
gularity at the center in the course of time. The central 
singularity is spacelike and not naked. In the case of flat 
or open universe models the singularity is weak and has 
distinct apparent and event horizons. The apparent hori- 
zons are not everywhere spacelike, to be compared with 
the Schwarzschild one which is null everywhere. This has 
already been noticed in a general context by [s^ . While 
the apparent horizon is defined by the surfaces R = 2M, 
similar to the Schwarzschild horizon, the even horizon 
is further away. Models we have proposed show that 
one has to expect new effects while considering dynami- 
cal cosmic black holes. The simple Schwarzschild static 
model may not reflect all the phenomena one may expect 
in observational cosmology, and the black hole thermo- 
dynamics. Even the simple concept of mass is not a triv- 
ial one in such a dynamical environment, the answer to 
these questions are beyond the scope of these paper and 
will be deal with in future publications. 



APPENDIX A 

The curvature function /(r) has a zero point where 
it changes sign for models III and IV, corresponding 
to two different solutions. Therefore, we have to take 
care of joining two solutions across the hypersurface 
defined by f{r) — to be continuous. This is done by 



looking at the metric functions and their derivatives to 
be continuous. 

Let us first look at the model IV. There we have to look 
at the metric function R and its derivatives, R, i?', R, R 
and i?', at the point r = 1 where / vanishes. From the 
following relations derived from the Einstein equations 

m 



R 



M 
R?' 



R 



_ M' MR' f 
RR RR^ 2i?' 



and 



M' 2MR' 



(Al) 



(A2) 



(A3) 



we infer that these second derivatives relevant for the 
Einstein equations to be continuous on the hypersurface 
/(r) = are continuous if the /, R, R' , R, M', and M 
are continuous. Now, because of the continuity of /, M', 
and M, we just have to prove the continuity of i?, i?, 
and R'. 

Let us look first at R and its derivative R' . In the case 
of r < 1 we have 



R = —^1 - COST]), 

1 — r 

(l-r)i-5 
b{r) 



(A4) 



where a{r) = + r^, h(r) = 1 + r^-^ , and a(l) — 2, 
a'(l) = 2.5, 6(1) = 2, 6'(1) = 1.5, and 



R = 



— r sinrj 

b 1 — COST] 



(A5) 



, _ a'(l -r) + a 



R' = 



(1 — COSlfj — 



(l-r)2 1-r 
sin'q 1.5(1 - r)°-56 + &'(1 - r)i-5^ 
1 — cosrj 

Defining 1 — r — x, we have 

77 - sim-j = — - 0{i]^) = —t. 



(A6) 



(A7) 



Therefore, to first order in 77 we have rj — \/ii^/x. Now 
taking the limit x — > 0^ we obtain 

2 

lim R{x) = lim — (1 — cos?7)= lim ( 0[ri^)/x) 



2/3 



(A8) 
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which is a well defined quantity. 

In the case of r > 1 we have 
a{r) 



and 



R 



r- 1 



{coshrj — 1) 



1)1.5 

smhrj — rj = — -r-r^ — 
6(r) 



a^/r — 1 sinhrj 
b coshrj — 1 ' 



(A9) 



(AlO) 



Therefore the continuity of R across r = 1 is established. 
Similar calculation for the first derivatives shows the con- 
tinuity of R' and R having well defined values on both 
sides of the r — hypersurface: 



R!{1) = 2.5(3i) 



2/3 



4(3^)1/3 



and 



(A14) 



R' = 



a'(r- 1) 



(r- 1)2 
sinhrj 1.5(r — 1)*''^6 



{coshrj — 1) 
b'{r 



coshrj — 1 
Now, defining r 



a 



1) 



r- 1 

1.5 
— t. 



52 

1 — and noting that 



(All) 



m = 



(3i)l/3 



(A15) 



sinhrj 



we obtain to first order of rj the relation rj — 
Therefore, 



lim R(x) 

X— 0+ 



lim — {coshrj 

r->o+ a; 



lim 



rr_ 

X 



(A12) 
"its/x. 



0{rj')/x 



= (3i)2/3 



(A13) 



The case of model III is similar except for the hyper- 
surface defined by g{r) = ^ = e"'' - ^^^20000 = 
the root of e"''" — 20000 ~ ^ being at a point r — ro 
different from r = 1. It is easy to see that g{r) is an 
analytic function at r = ro , and can be approximated by 
g{r) « g'(ro)(r — ro) -l- ^— j^(r — ro)2 + .... Similar calcu- 
lations verify the continuity of the metric function R and 
its relevant derivatives across the hypersurface r = ro. 
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